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We use non-linear gauge-invariant perturbation theory to study the interaction of an inflation pro- 
duced seed magnetic fleld with density and gravitational wave perturbations in an almost Friedmann- 
Lemaitre-Robertson- Walker (FLRW) spacetime with zero spatial curvature. We compare the effects 
of this coupling under the assumptions of poor conductivity, perfect conductivity and the case 
where the electric field is sourced via the coupling of velocity perturbations to the seed field in 
the ideal magnetohydrodynamic (MHD) regime, thus generalizing, improving on and correcting 
previous results. We solve our equations for long wavelength limits and numerically integrate the 
resulting equations to generate power spectra for the electromagnetic field variables, showing where 
the modes cross the horizon. We find that the interaction can seed Electric fields with non-zero curl 
and that the curl of the electric field dominates the power spectrum on small scales, in agreement 
with previous arguments. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 



fvq I Large scale magnetic fields of varying amplitudes are present in entire galaxy clusters, individual galaxies and high 

redshift condensations. Such fields are observed on characteristic scales of '--^ 1 Mpc and arc of micro-Gauss strength, 

'TT' 10^^ — 10^^ G [il, Hi- Despite their ubiquity, their origin is still a mystery. There are literally tens of candidate 

j^. mechanisms proposed to explain the origin and evolution of such fields, spanning different theories of physics [3|. It is 

I ■ now widely believed that the structure of magnetic fields in spiral galaxies is consistent with the dynamo amplification 

wp,, mechanism. The dynamo mechanism can produce amplification factors of up to ^ lO'^ but requires a seed field in 

i_i ■ order to operate and thus cannot explain the origin of magnetic fields. Additionally, adiabatic contraction of magnetic 

[ flux fines during structure formation can enhance galactic fields by a factor of ^ 10'^. 

\^ ' Among the physical mechanisms proposed to explain the origin of the seed field is one due to Harrison [J] . This 

J^ ' mechanism rests on the fact that non-zero vorticity in the pre-recombination photon-baryon plasma can generate 

jy-s . weak magnetic fields of about ~ 10^^^ G. However, vorticity is not a generated mode at first order in perturbation 

f~*) ■ theory and has to be put in as an initial condition. Second order treatments of the pre-recombination plasma in 

\^ , terms of a Kinetic theory description has also been used to generate the required seed fields [5l-[ig|. The key idea 

-y! ■ is a preferential Thompson scattering of photons off free electrons, over the scattering off protons (the scattering 

f^ ' off protons is suppressed by a factor {nie/mp)'^) which induces differences in the proton and electron velocity fields. 

^Sj ■ Electric fields are then induced to counter charge separation between the electrons and protons. The generated electric 

fields will then feed in the magnetic induction equation to generate magnetic fields at second order in perturbation 

theory. The photon anisotropic stress also couples to the electron velocities and contributes to the magnetic field 

sources. In addition, other arguments relying on electroweak phase transitions [ll|, [l2| , topological defects [1^ , velocity 

/\ . perturbations [IJ] etc. have been proposed as candidate mechanisms. The generated fields, however, are usually too 

H ' weak to leave any detectable imprint on the CMB Q . This is not surprising given the form of the fluid quantities of 

- - -' a magnetic field. In particular, the energy density hb = ^^/2, the isotropic pressure pb = B^/6 and the anisotropic 

pressure H^b = B/^Bi,\ of a field generated at second order will manifest at fourth order in perturbation theory, which 

is not relevant for CMB anisotropics. 

In addition to meeting the right strengths, the generated fields must be of the right scale to match those observed 
today. One of the problems of primordial generation mechanisms in general is that although some may reach the 
required strengths, they are causal in nature. This means that their coherence scales cannot exceed the Hubble scale 
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during the time of magnetic field generation. By comparison, the galactic scale today is well outside the Hubble 
scale at such early epochs. Moreover, the small scale fields i.c, those that are already sub-horizon before matter- 
radiation equality cannot reach the recombination epoch due to micro-physical mechanisms such as magnetic and 
photon diffusion processes 0]. 

Inflation and other pre-Big Bang models capable of causally producing super horizon perturbations are often 
invoked to circumvent this scale problem. However, the residual magnetic fields surviving the exponential expansion 
accompanying many inflationary models are thought to be too weak to be of cosmological relevance [65|. New physics 
often has to be introduced such as exotic couplings of the electromagnetic field to other fields such as the dilaton 
field to avoid the accompanying exponential dilution of the magnetic fields [la]. The primordial fields are also 
constrained by the fact that the anisotropic stress of the produced magnetic fields contains a spin-2 component and 
will result in an overproduction of gravitational waves at horizon crossing which is inconsistent with standard Big 
Bang Nucleosynthesis constraints [la, [I4I . 

Apart from studying the generation of magnetic fields, one can also study interactions of a pre-existing magnetic 
field with gravitational degrees of freedom. This is often studied in the context of amplification of the seed magnetic 
field or gravitational wave detection. Much progress has been made in this area [l8i - l2ll | . Most of these studies however 
have been restricted to focusing on the interaction of magnetic fields with tensor perturbations; In this work we revisit 
and extend the work presented in [la [20: 122]; to include scalar perturbations in the matter fluctuations. 

When using perturbation theory about a FLRW background to study the interaction, one is immediately faced 
with the problem of how to embed the seed magnetic field into the background. The isotropy of the FLRW spacetime 
does not readily allow for any direction preference that may be introduced by a vector field. There are several ways to 
handle this and we mention briefly just three of them. One can treat the seed magnetic field as a zeroth order quantity, 
subject to the assumption that the energy density of the field be small compared to the energy density of matter 
B^ <^ ^ and that the anisotropic stress is negligible liab ~ Bi^^B^,^ k, 0. With these approximations, the energy 
density of the magnetic field cannot alter the gravitational dynamics of the background spacetime; this approach 
is often referred to as the weak-field approximation. Another approach is to treat the seed field as a statistically 
homogeneous and isotropic random field with (B) = but {B"^) ^ and so, the seed field does not introduce any 
directional dependence in the background spacetime. One can then easily employ statistical methods to quantify 
the field's behavior. Another possibility is to leave the background spacetime untouched but treat the seed field 
as a first order perturbation, using a two parameter approximation scheme to characterize the perturbations in the 
electromagnetic and gravitational field; this is the approach we adopt in this work. 

One can go a long way in comparing the different perturbation schemes. For example, in the weak-field approxima- 
tion, the induced magnetic field will be at first order, a well understood regime in perturbation theory. While in the 
two parameter case, the induced field will be at second order [63, a regime that is not so well developed. Nevertheless, 
for the purposes of our work, the two approaches are mathematically equivalent. The apparent differences between 
them is as a result of relabeling of spacetimes, i.e. 'First order' in the weak-field approximation corresponds to 'second 
order' in the two parameter case. Indeed, Maxwell's equations and thus the Einstein-Maxwell system takes the same 
mathematical form in both of these approaches . They both use the machinery of relativistic perturbation theory and 
are thus prone to gauge issues, see |23|, 12J] for example. 

The present article is structured as follows: we present details of our pcrturbativc framework in § IIIIl After a 
presentation of the interaction equations in WI[ we present the derivation of the equations describing the induction 
of EM fields in ^3 VII Al and I VII Bl for a general current and a note on how to evaluate the induced electrical current in 
WII CI We present the power spectra of the induced magnetic field variable in ^X] and finally a summary in IXII We 
employ the 1-1-3 covariant approach to perturbation theory [25| and follow [26| by adopting the more geometrically 

motivated metric signature ( h ++) and we use geometrized units SttG = c = 1, where G is the gravitational 

constant and c is the speed of light in vacuum. 

II. PRELIMINARIES 
A. 1+3 spacetime splitting 

One of the nice aspects of the l-f3 covariant approach to General Relativity (GR) is that the underlying dynamical 
equations have a stronger appeal from a physical point of view, as compared to the quasi-linear, second-order partial 
differential equation form, which the EFE take in the metric based approach. 

The approach is based on a 1-1-3 decomposition of geometric quantities with respect to a fundamental four velocity 

dx"" 

Ua = -J-, UaU" = -1 , (1) 

(XT 



where x"" are general coordinates and r measures the proper time along the world line. The key equations governing 
the full structure of the spacetime are derived from the Ricci and the once and twice contracted Bianchi identities 
applied to the 4- velocity vector [23 ■ This splitting uniquely defines two projection tensors 

U\ = -u'^Ub => U\U\ - U\, U\ = 1, Uabu'' = Ua , (2) 

Kb = (Jab + UaUb => Kc^b = Kb, K^ = 3, habu'' = , (3) 

which project along and orthogonal to the 4- velocity w°. We define two projected covariant derivatives, the convective 
time derivative along u° and the spatially projected covariant derivative 

Q-V.., ^ w'^Veg^-V.., and B,Q'^-\...,^h%---h\h\---K^hf,yfQP-\...^ (4) 

respectively. The basic equations are then characterized by the irreducible parts of the first covariant derivative of Ua 

VaUb = -UaAb + T^aUb = -UaAb + T^^Kb + (Jab + ^ab , (5) 

where Ab = u'^'^aUb is the relativistic acceleration vector representing the effect of inertial forces on the fluid; 
DqM° = is the rate of volume expansion; aab = D{aWb) is the symmetric trace- free rate of shear tensor, describing 
the rate of distortion of the fluid flow; cuab ~ D[a'"fc] is the antisymmetric vorticity tensor, describing the rigid rotation 
of the fluid relative to a non-rotating frame. 

B. FLRW background 

We choose as our background the FLRW models, which are spatially homogeneous and isotropic. Thus, relative 
to the congruence Ua, the kinematical variables have to be locally isotropic, which implies the vanishing of the 4- 
acceleration Ua = 0, the rate of shear aab = and the vorticity vector Wa = 0. Spatial homogeneity implies that the 
spatial gradients of the energy density /x, pressure p, and the expansion Q vanish, i.e Da/i = D^p = D^O = 0. Moreover, 
the FLRW spacetime is characterized by a perfect fluid matter tensor, i.e tt ~ qa = 0. These restrictions imply that 
the spacetime is conformally flat, i.e the electric and magnetic parts of the Weyl tensor vanish, Eab = Hab = 0. This 
leads to the key background equations, the energy conservation equation 

/i= -(l + w)e/i, (6) 

the Raychaudhuri equation 
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e = --e^--M(i + 3w) + A, (7) 



where w = p/ fi and the Friedmann equation 



^, + A^Iq' + ^. (8) 



III. PERTURBATIVE FRAMEWORK 

As already mentioned, a FLRW spacetime cannot readily host magnetic fields, as their anisotropic stresses 
Hafo = B/^^Bb) 7^ will break the isotropy. We thus treat the background magnetic field Ba as a first order per- 
turbation to the isotropic spacetime. This lends the energy density, the isotropic and anisotropic pressure of the field 
to second order in perturbation theory. 

We then proceed by adopting a two parameter perturbative framework [28l432l |. Fundamentally, this consists of 
separately parametrizing the gravitational and Maxwell field perturbations in two expansion parameters €„ and es, 
representing the amplitudes of the gravitational and electromagnetic field perturbations, respectively [Wl . l22l . |2^ . 
Using this parametrization, any quantity Q" in the physical spacetime can be expanded in the form, 

(0,0) (1.0) (0,1) (1.1) 

Q-... = 44 Q-... +44 Q-... +e°4 g-... +44 g-... +o(4 , 4) , (9) 

where the first term on the right represents the background term; the first and second terms represent the first order 
gravitational and electromagnetic perturbations respectively; the fourth term represents the non-linear coupling we're 



looking to investigate; the higher order terms represent self-couphng terms of order e™ and e^, m,n > 2. In general, 
terms describing the coupling will be of the form e^e^, where, in this work, we restrict the perturbative order to 
C(eje^) and therefore neglect terms of order ©(e^e^), 0{e^e^) and higher, resulting from the self-coupling of the 
fields; this includes gravitational couplings with the magnetic anisotropy Xlaf, = —Bi^aBb), leading to ©(e^e^) terms. 
We will generally refer to quantities of order ©(eie^j) simply as non-linear and reserve the designation 'second order' 
for terms that are of order e^ and e^- . As in [l^, [22: [23 1 ^^^ '^^^ visualize this framework as a hierarchy of spacetimes 
to label the different perturbative orders. 

We make the common assumption in the literature that the perturbed spacetimes have the same manifold as the 
background spacetime i.e. we consider the perturbations as fields propagating on the background spacetime [321 l33l|. 
In this treatment, therefore, we restrict the possibility that the perturbations may alter the differential structure of 
the background manifold and so we neglect issues of backreaction. 

We're also interested in studying this coupling in a gauge-invariant manner. The gauge problem in relativistic 
perturbation theory has been dealt with in the literature, see for example [32, I34l-i37i . The Stewart & Walker Lemma 
[37| serves as a basis for the generalization of gauge invariance to arbitrary order [S^, |3g . It follows that a quantity Q 
is gauge invariant at order 0(e™e^) if and only if Q'"' and its perturbations of order lower than 0(e™e^) are either 
vanishing, or a constant scalar or a combination of Kronecker deltas with constant coefficients |28l . [32| . 

Since the interaction terms are of order C(e^eg) we have that the induced magnetic field Ba will be of the same 
order; we also assume that the electric field Ea will be of the same order as the induced magnetic field. Clearly Ba 
does not satisfy the criteria for gauge invariance at ©(eje^) since it is neither vanishing nor a constant scalar at 

©(e^e^). To this end, we make use of the same auxiliary variable Pa = Ba + |Si?a identified in [l^, [22, [23. We 
do not however integrate /3a to recover the gauge-dependent magnetic field, but treat it as the fundamental variable 
whose deviation from zero quantifies deviation from the adiabatic decay of the magnetic field. 

IV. THE EINSTEIN-MAXWELL SYSTEM 

The Einstein-Maxwell equations lAlOl contain terms that couple the electromagnetic fields to the gravitational fields. 
These can be written at ©(e^e^) by discarding higher order terms. This results in the two propagation equations, 

B(a) + \QBa = aabB'' - curl Ea , (10) 

E(a) + \QEa = curl Ba + EabcA'' B^ - Ja , (11) 

subject to the constraints, Da-E° = = Dai?". Following [33 we make the following comments: (i) The magnetic 
field Ba appearing in the equations [T3| and [T4| multiplied by the gravitational variables should not be the same as 
the Ba appearing alone. The variable B'^ is a mixture of linear and non-linear quantities (the seed magnetic field 
and the induced field) while the terms involving i3" are a product of first order quantities. One has to keep this in 
mind when integrating the equations, (ii) The system is not gauge-invariant as already mentioned in § IIIII This can 
be attributed to the mixture of linear and non-linear terms in the system. In the covariant approach to perturbation 
theory, the solution of perturbed differential operators is never sought, one can get around this by making sure that 
the differential operators involved operate on quantities of the corresponding perturbative order. 

In an attempt to cast it in a consistent and gauge invariant manner, we introduce the following non-linear variables: 
The fundamental variable f3a measuring deviation from adiabatic decay, la describing the interaction with shear 
distortions and ^a describing interaction with density perturbations. These are defined as, 

Pa = B^a)+^eBa, h ^ (JabB"" and ia ^ tabcA'' B" (12) 



and results in the following system. 



3 

where we have written curl Ea = Sa and curl Ba = ^a for brevity. 



Pa = la-Sa, (13) 

E(a) + ^QEa = ^a+^a-Ja, (W) 



V. THE LINEAR EQUATIONS 

A. The linear magnetic field: 0{eg) 

We treat the seed magnetic field as a first order perturbation to the spacetime. The seed field may have its origins in 
inflation or other mechanisms based on string cosmology, in which electromagnetic vacuum fluctuations are amplified 
due to a dynamical dilaton or an infiaton field [la]. We assume that at order ©(e^e^) the electric fields are small 
compared to the magnetic fields, i.c £'^ <C -B^. Thus, in the absence of diffusive losses or amplification, the induction 
equation 1131 takes the frozen-in form, 

B(a) + \<dBa = , (15) 

regardless of the equation of state or plasma properties of the cosmic fiuid. It follows then that the magnetic field 
decays adiabatically as Ba oc a~^, where a is the cosmological scale factor. This adiabatic decay arises from the 
expansion of the Universe which conformally dilutes the field lines due to flux conservation. The frozen-in condition 
[T5]does not discriminate between homogeneous (DaS^ = 0) and inhomogeneous {T>aBb ^ 0) magnetic fields. For an 
inhomogeneous field the spatial gradients of the seed magnetic field T)i,Ba are of the same order as Ba and evolve as 
DbSa ex a"3. 

B. Gravitational perturbations: 0{eg) 

The Weyl tensor Cabcd represents the free gravitational field, enabling gravitational action at a distance. In analogy 
with splitting the Maxwell field tensor Fab vaio a magnetic and an electric field, Cabcd can be split covariantly 
into a 'magnetic' part Hab = ^^adeC'^'^fj^u'^ and an 'electric' part Eab = Cabcdu'^u'^- The electric part of the Weyl 
tensor describes tidal effects, akin to the tidal tensor associated with the Newtonian potential, while the magnetic 
part describes the propagation of gravitational radiation. The Weyl tensor vanishes in the conformally flat FLRW 
spacetime and so Eab and Hab are covariant first order gauge invariant (FOGI) quantities in the Weyl curvature. We 
also define the FOGI variables Xa = aDa/x and Z = aDaO to characterize density perturbations. Now, the system 
governing gravitational perturbations is given by the following propagation equations [67[ , 

2 

<^{ab) + o©^at = D(aA> - Eab , (16) 

Hi^ab) + QHab =- -curl Eab , (17) 

E(^ab) + 'dEab = curl Hab - -^/^'(l + W)<^ab , (18) 

X,^a)-QwXa = -{l+w)Za, (19) 

^"^3 2^ 3(l + w) V 3 ^^^ y "^ 1 + u; 

In addition to the propagation equations above, the following constraints have to be satisfied, 

2 1 

aD'^CTbc = -Zb, aD^Ebc = xA'-^b and Hab = curl CTah , (21) 

where we have set the vorticity to zero (wq = 0), see also [38|. Note that at first order in gravitational perturbations, 
the only source of vector modes is the vorticity uja', since, we neglect the effects of vorticity, Wa = 0, all the vector 
modes vanish. The shear tensor (Tab can then be irreducibly split into scalar and tensor contributions as [39| 

<^ab = cTab + <^lb where curl Uab = 0, and D^crjf, = . (22) 

The pure tensor modes can be used to characterize gravitational waves [40] • The scalar part of the shear couples to 
density perturbations and is related to the clumping of matter via the constraints 1211 

By differentiating [TBI and using [T51 and one of the constraints W\\ to substitute for Eab and Hab, one arrives at a 
forced wave equation for the shear, 



Z>a) + -QZa = ~-ilXa ~ ,,^ , , --6^ + ^l + K] Xa - -—B'^Xa . (20) 



2 5 . 



1^2 1 3 5, 

-9^ + -n »+ -A 

9 6^ 2^ 3 



<Jab = 



a2(l + w) 



Xab + xQ-^ab 



(23) 



where Xab = —(1 + w)a^D(a^(,)/w^ = oD{a'^b)- We need an evolution equation for Xab in order to close equation!^ 
One can start from [TO] and [5D] to write a wave equation for Xa then taking the comoving spatial gradient of the 
resulting wave equation will yield the following, 



Xn. 



wWX„ 



2\ • 1 

3 j e-Yab + -^(3u; + \)(w - \)Xab - 2u;AA'at = 



(24) 



In including scalar perturbations, we have explicitly coupled the shear tensor to density perturbations. This shows 
that density gradients source distortions in the Weyl curvature and vice versa. Hence, knowing the shear allows one 
to compute density gradients and knowing density gradients one can compute the scalar part of the shear [4l| . 



VI. THE INTERACTION: ©(e^e^ 



The Maxwell fields couple to Weyl curvature through the shear term and density perturbations through the accel- 
eration terms and the non-linear identity 1281 In the case of pure tensor modes in the shear tensor, the interaction 
variable la = <y^ab" ^^^ shown to satisfy a closed wave equation, for both a homogeneous |18| and an inhomogeneous 
[22l | seed field Ba- Here, wc include contributions from scalar perturbations in the shear, which give rise to source 
terms due to coupling with density perturbations. In this case la satisfies a forced wave equation, 



where the forcing term C^ is given by. 



13^2 1 5 7/ 

— 9 ii uw + -A 

9 6^ 2^ 3 



^0. — ^a 1 



[S^a) + e^a) 






To close the above system, we give the companion wave equation for Sa — aS''D(a<^6) ^s, 



Sa - wB^Sa + (2 - w)eSa 



^(1 - w){A + e^) + ln{l + 3w){3w - 5) 
3 6 



Sa=0. 



(25) 



(26) 



(27) 



We note, for later convenience (^ IVIIB[) that the forcing term C^ = in a matter dominated universe {w — 0) i.e la 
decouples from Sa when w = 0. 



VII. INDUCTION OF EM FIELDS 

We introduce non-linear gravitationally induced 'effective current' terms C^, Cf and C^ which are made up of the 
coupling between density and gravitational wave perturbations; these will act as driving forces of the induced Maxwell 
fields. 



The Electric field 



We show how the coupling of gravitational perturbations with the seed magnetic field can induce Electric fields. 
Here we give wave equations for the induced Electric field Ea and its rotation (oa. In deriving the wave equation for 
Ea, we differentiate [14] and equate the result to the non-linear identity. 



(curl Ba) = curl /3„ - 9 curl Ba + HabB^ 



1 



3a(l + u;) 



EabcB^ SwX" - 2X' 



(28) 



obtained from the commutation relations (Appendix lASp and we have used Equation [T2] to rewrite the acceleration 
terms. The resulting wave equation is found to be. 



El^a) - B^Ea + ^QEa 



^e' + ^fi{l-3w) + ^A 



Ea - Ca , 



(29) 



where C^ is a gravitationally induced source term given by, 



Cf = curl la + HabB'' 



1 



a(l + w) 



^abc 



w-'^yB'x^')- + eL-^]B'x^ 



^U a Ua ■) 



(30) 



and J7a is the 3-current. The terms involving eabc in C^ vanish when the magnetic field i?° is parallel to the fractional 
density gradient A"". Taking the curl of [IHl results in the equation governing the rotation of Ea, 



2 7 • 

©a ^ D ©a + —QS'a 



Ie» + iMi-9u.) + 5A 



€.. 






(31) 



where the source term Cf = curl C^ is given by, 



-(curl Ja)-- -ecurl Ja + 2D^D[,4] + eacdBbB' H'"' 



(32) 



B. The Magnetic field 

As already mentioned, the induced magnetic field will be characterized via the variable /3a = Ba + ^QBa- On using 
[T51 E5] and [211 one can write a second-order equation governing the evolution of the fundamental variable /3a ■ This 
can be written in either of two forms: in terms of la or (oa, corresponding to using [T3] as a constraint to either of 1311 
or [25] respectively. Recall that both la and S'a satisfy wave equations of the form C[Ia] = C^ and £[(o'a] = Cf , where 
the C*s are source terms. 

Using covariant harmonics [4l|, one can already notice from l25l and l31l that the eigenfunctions used to harmonically 
decompose la and S'a are not the same for a general perturbation [68| . Consider the induction equation 1131 and write 
it as Pa ~ 'Ylik^o-I(k) ~ Qa'o(fc)), whcrc Va and Qa are distinct eigenfunctions of the Laplace-Beltrami operator, i.e 
T^a 7^ Qa- For the separation of variables technique to work for Pa, one must eliminate either la = Val(k)i along with 
its source terms C^ or S'a = Qa^{k) along with its source terms Cf^ . In this way, Pa can then be expanded in terms of 
one set of complete eigenfunctions. This presents a problem: since both la and Sa are coupled to source terms C^ and 
Ca respectively at second-order, both C^ and C^ will still couple to the P a eq uation at this order, thereby introducing 
the differing set of eigenfunctions Va and Qa- A similar problem arose in |42| . due to the inclusion of a vorticity term. 

It is possible to do away with C^ in equation [25] by requiring that w = and this alleviates the problem [62| • We 
shall then henceforth restrict to the pressureless dust [w = 0) case and write the Pa wave equation in terms of Sa- 



where. 





ha) - D'/3a + 3e/3(a) + 


cf = 


^^64 + 


2^2 1 2/ 
-3^+3^-3^ 


Sa + (Cl 


- 


2 


-^(aSiaD^A-fe])-- ^e(aB[aD^Afb]) 



13^2 1 7/ 



Pa 



c. 



p 



(33) 



Sa + (curl Ja)- + -e curl Ja - 2 D^Dla/fc] - EacdBbB'H' 



db 



(34) 



Note that while we keep Sa ~ aB'^D/^a^t) distinct from aBjaD Xi,] in real space, their evolution equations can be 
made equivalent in harmonic space by a suitable choice of eigenfunctions [70| . We shall thus write S'(^) in place of 
B(n)X{k) to avoid introducing another letter to denote the latter. This should not lead to any ambiguities. 

C. Tiie Electric Current 

1. Limiting cases: poor and perfect conductivity 

To close the above system, one needs to take care of the current term Ja appearing in [301 [321 and [Ml This term 
depends on the electrical properties of the medium. It is given in terms of the Electric field Ea via Ohm's law, 



Ja — ^Ea, 



(35) 



where <^ is the electrical conductivity of the medium. In this section, we consider only the limiting cases of very high 
(^ — )■ oo) and very poor conductivity (<; — )■ 0). Under the assumption of poor conductivity, the currents vanish Ja = 0, 
despite the presence of a non-zero electric field. In this case, one solves equations [29l [31] and |33l with the current 
terms set to zero. At the opposite end, the case of perfect conductivity, the electric fields vanish and the currents 
keep the magnetic field frozen in with the fluid. In this case, the current term satisfies. 



(curl J,)-+ ^ecurl J, = 2D^D[„/fa] + (acdBbB^H^'' + ^ 



'^{aB[aB''Xt])' - leiaB^aB'Xi,^) 



(36) 



and 1291 and 1311 are no longer relevant. One can verify that using this relation reduces equation 1331 to /3a = la, as can 
be confirmed also from the induction equation [T31 

One can also invoke the magnetohydrodynamic MHD ap pro ximation, which is valid for cold plasmas (pressureless 
dust can be well approximated by a cold plasma treatment) 14411 . Co ld p lasmas have components with non-rclativistic 
velocities and are thus mathematically easier to deal with [22, |4a, 13 • We consider a two component electron- ion 
plasma and assume that the motion properties of the plasma on macroscopic scales are captured by the center of 
mass 3- velocity v'^ of the system i.e the difference in mean velocities of the individual species is small compared with 
the fluid velocity. We also assume charge neutrality of the cosmic plasma, i.e., the number densities of the electrons 
and ions rig and Ui are roughly equal, rie ~ nf, this guarantees the vanishing of the total charge pc = — e(ne — n^) « 
and the background 3-currcnt JJa ~ 0. In this case, the generalized Ohm's law is given by 



J{a) = <Ea + eab^B^) , 



Mefe 



fJ-iV, 



A^e +Mj 



(37) 



where the subscripts e and i denote quantities for electrons and ions respectively. The center of mass 3- velocity Va of 
the electron-ion plasma can be shown to satisfy the linearized Euler equation. 



V{a) + ■^'dVa = 



(38) 



In the ideal-MHD environment, the conductivity of the medium is very high (<j — > oo), then Ea + CabcV^B'^ — > in 
order to keep the current J'a finite. This readily gives the Electric field Ea and its rotation S'a as, Ea = —tabcV^B'^ 
and Sa = 2i3[aD v^. Using [TSl and [551 one can show that. 



E 



(a) 



QEa = 0, 



and 



f&n. 



-e^fa^o. 



(39) 



With these, the 3-current J'a satisfies. 



(curl Ja)-+ -Ocurl Ja = 2D%aIb] + eacdBbD'H" 



-'^{aB[aB'Xb])'- ^eiaB[aT)'Xi,f, 



(40) 



Substituting [40] into [33] results in 

/3(a) 



B^l3a 



ie/3ia) 



e' 



6 



M 



-A 



Pa 



B^£a 



--&' 



7 A* 



A 



<&n. 



(41) 



The application of the ideal MHD approximation in cosmology has often been criticized as rather being of practical 
appeal rather than of physical one [43| . Ideally, the curl of Ea should be the outcome of a rigorous treatment of the 
physics of the particle interactions in terms of a kinetic theory description, see for example [3, [To] ■ 

2. Intermediate case: Finite conductivity 

The case of poor conductivity may not be much relevant in the post recombination epoch as the universe then 
acquires very high conductivity. The perfect conductivity case, while relevant, may be thought of as an idealized 
notion. We thus turn to the finite conductivity case. The conductivity of the post decoupling era can be modelled by. 



<; = 



10"s-i 



mf,n^crT 



(42) 



where, Ue is the density of free electrons, e is the electric charge of an electron, rrie is the mass of an electron, n^ is 
the density of photons and ctt is the collision crossection. For a perfect fluid, the ratio n^/ue is constant, see [4^ for 
example. 

Assuming that Ohm's law holds (Equation I35p. we may write the current terms in 1341 as. 



(curl Ja)- + -ecurl Ja = '^^ + 3 O <r^a 



(43) 



where we have assumed that spatial gradients of the conductivity may be neglected {Da'i ~ 0) and that the conductivity 
is constant in time (■^ ~ 0). Substituting 1431 in the wave equation 1331 for /3a results in. 



/3(,) - B'^Pa + 3ei3(^a) + 



— 8^ 11 + -A 

9 6^ 3 



Pa = Cf 



(44) 



where the source term Cf is now given by, 



e 



e^a 



e 



u|e^4. 



-A 



f&n 



2D''D[,4]-e„,rfBbD^iJ'^^ 



|(a%D''Afb])-- ^e(aB[„D''A'b]) 



(45) 



Note that the electric currents J'a, electric fields Ea and the conductivity <^ are all simultaneously finite. The simplifi- 
cations that arise due to the characterisation of the limitting cases (j7a = for poor conducting mediums and Ea = 
for perfect conducting mediums) are no longer applicable in the case of finite conductivity. One then needs a proper 
model for the electric currents to ensure that the initial conditions for J7a and Ea are not chosen independently. There 
are several ways in which one can model electric currents, all resulting in terms of perturbative order el, see [8| for 
example. While these terms can be seamlessly accommodated in our framework, they have the undesirable effect of 
seeding magnetic fields. This will lead us away from the isolated effects of the amplification of an already existing 
field. Inclusion of such terms will therefore lead us to overestimate the effect of the amplification. With this in mind, 
we restrict to the limitting cases of lVIICTl 

VIII. THE INDUCED FIELDS 

We now treat separately the induction of electromagnetic fields due to interaction with scalar and tensor pertur- 
bations. To this end, we expand the perturbation variables in terms of a helicity basis (Appendix I A II) . In addition, 
we use a unified time variable whose defining equation is f = ^Hi instead of proper time, to re-write the relevant 
equations [7l|. We have to substitute for fx, 9 and a, appearing in the perturbed equations, from the zeroth order 
equations. We restrict our treatment to zero cosmological constant A = and flat spatial sections K = 0. Friedmann 
equation then reduces to /z = 6-^/3, where 6 is given by O = SHi/r; the scale factor a evolves as a = air"^/^. 



A. EM induction due to scalar perturbations 



In this case, the coupling of a seed field with gravitational perturbations is described by the variable la and Sa] 
these variables become sources of electromagnetic fields. 



• Interaction terms: Equations [^ and [?7l for the interaction variables la and Sa, respectively become, 

:/".> + =-I'^^ + ^I(e^ = 0, (46a) 



9 27 25 



997 



'(£) 







(46b) 



Note that since w = 0, the entire system has decoupled from aB'^T) /a'^b) : however we still need an equation 
for S(^i) because of the coupling with ai3[aD X^,^ in Equations [32] and [34l These interaction variables have the 
general solutions. 



/m(T)=Cir"i"/3 + C2r-5 



/3 



and 



'W 



iCgr-^/^ + ic4r-2/3 



(47) 



where the C^'s are integration constants. 



• EM fields: Equation [551 for the Electric field Ea becomes, 






aiH, 



.-4/3 



(fc^hn). ^-2/3^ 1 o/ ^ 



^(^'~"^M^'(^) 



i/. (^'^3il,r 



5, 



(«) 



It is much easier to solve for /3(£) from the induction equation 



ki)T 



273 -^W 



10 



(48) 



(49) 



once III) and Ei£\ are known, rather than from the wave equation [ 



B. EM induction due to tensor perturbations 



In this case, the transverse and trace-free parts of the shear tensor Uab characterize gravitational waves. The 
interaction with a seed field is then purely described by the variable /„ without any contribution from cither density 
or velocity perturbations. The generalized Ohm's law |37] in the MHD approximation also reduces to [35] We thus 
only need the equations for /3a, la and Ea- 



• Interaction variable: Equation [^Sj for the interaction variable la becomes, 

"2 



Q 97 

-T" A-—T' 4- 

4 W + 2r W + 



aiH, 



-4/3 



with the general solution. 



-5/2 



CiJi 



25 
2^2 



C2J2 



'(£) 



£ 



where Ci and C2 are integration constants, Ji and J2 are Bessel functions of the second kind. 
• EM Fields: Equation [25|for the electric field variable Ea becomes, 



(50) 



(51) 



4 (^) 27 ('') 



and we once again determine /3^i) from 



instead of using the wave equation (331 



i 



aiHi 



.-4/3 



E^e)=± 



C^k + n) 2/3 



hi) T rT57^^(«) 



(52) 



(53) 



IX. INITIAL CONDITIONS 

We need initial conditions in order to solve the equations in the previous section. The conditions are adapted as 
follows: for Pa we invoke Maxwell's equation [T^ 



Pa — la ^ 'S'a, Pa — la ^ <^a 

For the interaction variable la , we use the definition [T^] and Equation [T5| 

In = (JnhB In = (TnhB + (TnhB Bn = 



-QBa 



(54) 



(55) 



For the rotation of the Electric field S'a [72|, we use Maxwell's equation [T4l and the commutation relation [A13I to get, 

4. = -e^a+UabB^'-D^Ba (56) 
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where in this case i?a (without the tilde) is the induced magnetic field, and we have written the first order perturbed 
3-ricci tensor TZab as [2^ [23 



T^ab — ~^^(ab) ~ ^O^ab 



(57) 



We require that the gravitationally induced field variables Ea (and hence (oa) and Ba be zero initially. This leads to 
the following initial conditions for the perturbation variables: 



/: 



W 



/3, 



"U^U 



w 



/' 



w 



j(f)-°(n) 3"(fc)-°(n) 



'^w - 



<C', 



w 



'^(fe)^(n) 



(3' 



m 



-2S, 



w 



(< 



w 



2o-.i(fc)-S(„) + l<^\k)B\^) 



^-U BU 



(58) 



2S, 



W 



Following [2l|, |43, |43, we adopt the initial condition for the shear from {a/H)i ~ 10 ^. We choose the seed field to 
be B^ = 10"^" G, as typical of those produced around the recombination era [lOl |. 



X. RESULTS 

Given the system of initial conditions 1581 one can notice that the interaction variable la plays the fundamental role 
in the interaction process. If we set /q = initially, then no amplification takes place. We show the time evolution 
Ia{T) in Figure [1] on a log-log scale. A noteworthy feature is the rapid decay of la for both scalars and tensors. 
Although the interaction with scalar perturbations decays slightly slower, it essentially follows the same trend as the 
interaction with gravitational waves. We are thus led to conclude that even including scalar perturbations in the 
interaction, we reach the same conclusion as [l9[ and |45| that there is no significant amplification of electromagnetic 
fields coming from the interaction. 




FIG. 1: Time evolution of the interaction variable in log-log axes, 
slower than for tensors. 



Note that for the interaction with scalars, the decay is slightly 



The effect of the gravitational perturbations on the interaction is thought to be largest at the point where the 

2(b)| 



modes enter the horizon. This is clearly evident in Figures 2(a) and 
from Figure 2(a) One is that the spectrum for the interaction varia 



A couple of features are worth noting 
ible mimics that of gravitational waves. It is 
also consistent with the fact that gravitational waves start oscillating at horizon crossing. This is to be expected 
since although for a spatially inhomogeneous magnetic field Ba, the product /(£) == S(n)f(A:) becomes a convolution 
in Fourier space, I{k) ~ J^n B{n)a{k — n), we have only considered the mode-mode coupling case, I{k) = B{k)a{k). 
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K^ 




10""° V 
10 



10^ 



10^ 



10^ 



10' 



(a) 



(b) 



FIG. 2: Plots of Power vs scale (£); we define the power as Px = |a;(^)'^|. (a): Power spectra of the magnetic field variable I3(^i-j 
(green, solid), and the interaction variable I^i) (blue, dashed) at redshift z — for the tensor case, (b): Power spectra of the 
magnetic field variable /3(f) (green, solid), and the interaction variable I^i) (blue, dashed) at redshift z = for the scalar case. 

The power spectra for the case of interaction with scalars are not as interesting. There is no scale dependence on 
the interaction variable la , cf equation I46al This is because the Laplacian term for scalar perturbations conies from 
the acceleration vector which is identically zero in the dust case A = 0. 

It would be interesting to generalize our treatment to include the case of non-zero pressure. This will lend us to 
the radiation dominated era where one can incorporate photons in the plasma and can consider collisional effects as 
was done in [3, \l^ for example. One could treat the interesting case of simultaneous generation and amplification of 
magnetic fields by coalescing these phenomenon. 



XI. CONCLUSIONS 



We have carried out an analysis of the coupling between gravitational perturbations with electromagnetic fields 
as a possible means for magnetic field amplification. This carries to completion the work began in [1^ and [221. I^^ 
agreement with the work of [l9[ and [45| we argue that there is no significant amplification resulting from the interaction 
of magnetic field with gravitational waves. Even with the inclusion of density perturbations, the induced fields may 
still be orders of magnitude smaller. This justifies the perturbative treatment and our neglect of backrcaction. 

The induction of electromagnetic fields due to the interaction of a test magnetic field with gravitational waves was 
studied in [20[ using the weak-field approximation. We included this study here treating the background magnetic 
field as a first order perturbation and recovered similar results. This shows that there is no fundamental difference 
between the two approaches, apart from a labeling of spacetimes, which should not affect physical results. We also 
extended this study by using a proper non- linear perturbative framework. This framework was applied in [l8|, but 
an erroneous argument there led to the neglect of the rotation of the Electric field, thus restricting the study to 
perfectly conducting environments. This was refuted in j52| . In fact, upon inspection of[SS]one can conclude that 
even if one initially sets the rotation of the Electric field to zero, S'i — there arc non-zero terms on the right hand 
side of the initial conditions for S" that will seed a non-zero S". Wc also carry to completion the work in [22[ by doing 
a proper extraction of the scalar and tensor modes and numerical inte grat ions. In terms of the conductivity of the 
cosmic medium, |20[ restricted their study to poor conducting mediums, |18| to perfectly conducting mediums and [22| 
treated the MHD approximation. We carried our analysis for all three cases. We find that for tensor perturbations, 
the ideal MHD approximation is just the same as the perfect conductivity assumption of the fluid treatment. For 
scalar perturbations, we find an additional source term in the induced field (compared with perfectly conducting 
environments) due to the coupling of the seed field with scalar velocity perturbations. The current term J7a was 
neglected at all orders in |18j . in an attempt presumably to uphold the background magnetic field's homogeneity 
condition DaBt ~ 0. However, this is not necessary since introducing the current term at the non-linear order does 
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not break the condition DaBi, = 0. Also, one cannot consistently invoke Ohm's law for poor and perfect conducting 
environments without a current term. In |22| . an inhomogeneous seed field was assumed thereby requiring a first 
order current J'^ = p^v^ + pivf = —e{neV^ — Uivf) to uphold the condition Da-Bb 7^ 0. However, after decoupling 
(which is the era considered there), Thompson scattering is no longer efficient. Thus electrons and ions are tightly 
coupled by coulomb scattering at first order. Their velocity fields are therefore equal as they form a perfectly coupled 
baryon fluid [50, l5l[ ■ There can be therefore no currents generated at this order and the condition curl Ba ~ J a will 
render the seed field homogeneous. 

Both [l^ and [23j integrate /?a to recover the amplified magnetic field, after specifying a frame m°. While this 
takes into account the frame dependence of the magnetic field Ba-, it invalidates gauge invariance as the recovered Ba 
remains gauge dependent and takes the same value and form as it would have without the introduction of /3a- This 
is already pointed out in [52|, See also [5j] and [5J|. We do not integrate /3a but simply note that one can assign a 
physical meaning to the magnetic field variable /3a by noting that /3a = describes the background adiabatic decay 
of the fields. Any deviation from /3a = would then imply amplification of the background field. Moreover, /3a is 
a linear combination of terms that source magnetic fields through the induction equation [T31 Thus we can estimate 
the relative importance of each source term through /3a without having to integrate it to recover the gauge-dependent 



_B°. For example, we see from Figure 2(a) that the rotation of the electric field dominates at small scales compared 
to the interaction term. Observations of cosmological magnetic fields are difficult enough as it is, a new cosmological 
observable would lead to better understanding of studies in magnetic fields. While /3a may not be that quantity, it 
does arise naturally from Maxwell's equations. 

Also, one can readily write our key equations in terms of metric variables by adoption of a suitable tetrad as was 
done in [5[. 

Mechanisms that seek to generate magnetic fields, relying on non-linear perturbation theory are attractive for 
several reasons [53 • Among these is that they can easily blend in with known physics as they become relevant around 
the recombination era. This makes it possible to quantitatively evaluate the generated fields using CMB constraints. 
Progress in non-linear perturbation theory will allow us to investigate these non-linear effects in a manner that is free 
of spurious gauge modes [38|, |56[ . 
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Appendix A: Appendix 

1. Harmonic Splitting 

It is standard to decompose the perturbed variables harmonically in Fourier space; separating out the time and 
space variations [4l|, [5^, [53 • The idea is to expand the quantities in terms of eigenfunctions of the Laplace-Beltrami 
operator. To this end, we introduce the helicity basis vectors e'"\e<°' and e<+' defined by 

ei*'=-^(ei±*e^) (Al) 

where (e-'-,e^, k) form a right-handed orthornomal system with 62 = k x ei and we align e° with k. 
Using this basis, the scalar harmonic functions are given by, 

g(0) ^ ^k,xi (^2) 

Scalar type components of vectors and tensors are expanded in terms of harmonic functions defined from Q^^' as 
follows, 

QT = -7DaO(")-azfcae''=^-\ (A3) 

k 

Q^al = ^D{aD,)Q(°) = -a^ (kak - ^a,^ 6^"^^' . (A4) 
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Vector harmonics are given by 

Qi"' = ei±'Q(°) , (A5) 

Qit' = -^D(aelf 0^°) = azfc(,e<f e^^^-^' . (A6) 

While tensor harmonics are defined as, 

Qtl = vl^i*'er'Q(°) . (A7) 

2. Mcixwell's Equations 

The Maxwell field tensor Fab decomposes relative to the fundamental observer as, 

Fab = 2u^aEb] + eabcB" , (A8) 

where Ea = FabU^ and Ba = \e.abcF^'' are respectively the Electric and Magnetic field as measured by the fundamental 
observer moving with 4-velocity m". These are 3-vectors on the spacelike hypersurface, Eau"" = = Bau""- The 
Maxwell's equations are given by 

V[,Fbe] = and V''F,b = Ja , (A9) 

where J is the 4-current. These equations can be decomposed covariantly into the following |25ll60ll6l| 

ij^a) - curl Ba = -^eEa + aabE'' + EabM^ B'' + io^ E') - ^tQ ^{a) , (AlOa) 

B^a) + curl Fa = -l^Ba + (TabB^' + EabciA^ E' + Uj"" B') , (AlOb) 

= DaE" - 2ujaB'' - — , (AlOc) 

eo 

= DaB"" + iuaE" . (AlOd) 

The EM fields are solenoidal in the absence of gravitational vector perturbations. 

3. Commutation Relations 



(D„/)-i = D„/-ieDa/ + /A, (All) 

{'DaVb)l_ = ^aVb - ^Q^aVb - <D,Fb + ebcdHa'V + AaVb , (A12) 

(curl K)± = curl K - ^6 curl K - Cabc^'^'Drfl/^ + EabV'' - ^^eabcV'' A^ , (A13) 

curl curl 5„6 = -D^^^,, + (m + A - ie2)5a6 + ^D^^D-^^,), . (A14) 
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